Bose-Einstein condensates with attractive 1/r interaction: The case of self-trapping 
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Amplifying on a proposal by O'Dell et al. for the realization of Bose-Einstein condensates of 
neutral atoms with attractive 1/r interaction, we point out that the instance of self-trapping of the 
condensate, without external trap potential, is physically best understood by introducing appro- 
priate "atomic" units. This reveals a remarkable scaling property: the physics of the condensate 
depends only on the two parameters N 2 a/a u and 7/iV 2 , where N is the particle number, a the 
scattering length, a u the "Bohr" radius and 7 the trap frequency in atomic units. We calculate 
accurate numerical results for self-trapping wave functions and potentials, for energies, sizes and 
peak densities, and compare with previous variational results. As a novel feature we point out 
the existence of a second solution of the extended Gross-Pitaevskii equation for negative scattering 
lengths, with and without trapping potential, which is born together with the ground state in a 
tangent bifurcation. This indicates the existence of an unstable collectively excited state of the 
condensate for negative scattering lengths. 

PACS numbers: 03.75.Hh, 34.20.Cf, 34.80.Qb, 04.40.-b 



I. INTRODUCTION 

Bose-Einstein condensation of dipolar gases has at- 
tracted much attention in recent years P, S, H, 13, Q 
because it offers the opportunity to create degenerate 
quantum gases with adjustable long- and short-range in- 
teractions, which gives rise to a wealth of novel phe- 
nomena [y, 0, H, [s[. In particular, the achievement of 
Bose-Einstein condensation in a gas of chromium atoms 
(To| . with a large dipole moment, has opened the way to 
promising experiments on dipolar quantum gases [111 ]. 

As an alternative system with tunable interactions, the 
Bose-Einstein condensation of neutral atoms with elec- 
tromagnetically induced attractive 1/r interaction has 
been proposed. Here a monopolar, "gravity-like", long- 
range interaction, in addition to the short-range (van- 
der-Waals-like) interactions, takes the place of the dipole- 
dipole interaction in dipolar gases. A monopolar quan- 
tum gas could be realized according to O'Dell et al. [IH 
by a combination of 6 appropriately arranged "triads" 
of intense off-resonant laser beams. In that arrange- 
ment, the 1/r 3 interactions of the retarded dipole-dipole 
interaction of neutral atoms in the presence of intense 
electromagnetic radiation are averaged out in the near- 
zone limit [H, [lH, while the weaker 1/r interaction is 
retained. The resul ting atom-atom interaction potential 
in the near-zone is (PSl 



V u (r,7) 



with u = 



11 Ik 2 a 2 



4-7T cei 



Here, a(k) is the isotropic, dynamic, polarizability of 
the atoms at frequency ck, and / the intensity of the 
radiation. The quantity u determines the strength of 
the "gravity- like" interaction. The estimate for u given 
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by O'Dell et al. [12( for CO2 laser light of intensity 
/ = 10 8 W/cm 2 is equivalent to the attraction of two 
opposite equal charges with q w e/2000. However, by 
contrast with the van der Waals interaction, the 1/r po- 
tential acts over the entire sample, and therefore its con- 
tribution to the energy can become important. Instead 
of 6 triads of lasers a different arrangement with three 
rotating lasers has been proposed [l5| . 

Even though the experimental realization of such con- 
figurations is not yet at hand, the theoretical issues as- 
sociated with monopolar degenerate quantum gases are 
worthwhile investigating. In particular, as pointed out 
by O'Dell et al. [T4I, the intriguing new physical fea- 
ture that emerges is the possibility of self-trapping of the 
condensate, without external trap. 

In the theory of trapped Bose-Einstein condensates it 
is common to introduce as natural units for energy and 
length the quantum energy Hujq and the oscillator length 
ao = yJh/mu)Q of the trap potential. In the case of 
self-trapping, however, where the trapping potential is 
switched off, Hujq — > and ao — > 00. Thus these quan- 
tities become "bad" units. As a consequence, in their 
study of the physical conditions necessary to observe the 
transition from external binding to self-binding Giovan- 
nazzi et al. [l6| used the laser wavelength and energy as 
units of length and energy. 

It is the purpose of this paper to re-analyze Bose con- 
densates with attractive 1/r interaction using appropri- 
ate "atomic" units. This will first reveal remarkable scal- 
ing properties of the condensates. Next we solve the 
extended Gross-Pitaevskii equation for monopolar quan- 
tum gases numerically and compare with previous vari- 
ational results. Last, as a novel feature, we point out 
that our numerical calculations reveal the existence of a 
second solution of the Gross-Pitaevskii equation for neg- 
ative scattering lengths, which is born together with the 
ground state in a bifurcation " out of nowhere" . The ex- 
istence of the second solution indicates the existence of 
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an unstable collectively excited state of such condensates 
at negative scattering lengths. 



II. NATURAL UNITS, SCALING PROPERTIES 
A. The general case 

We analyze the physics of trapped monopolar gases, 
and in particular the limit luq — > 0, in terms of natural 
"atomic" units. From the analogy u e 2 /Atteq we can 
define a "fine-structure constant" 



u/hc 



(1) 



and can construct a " Bohr radius" and " Rydberg energy" 
in the usual way from the Compton wavelength Ac = 
h/mc and the rest energy mc 2 via 



Ac 



mu 



E„ 



2 2 

at mc 



2ma 2 



(2) 



Measuring lengths in a u and energies in E u , we can write 
the Hartree equation of the ground state of a system 
of N identical bosons in an isotropic external trapping 
potential Vo(r) = mui 2 r 2 /2, all in the same single-particle 
orbital ip, interacting via V u and the s-wave scattering 
pscudopotential V s — AiraH 2 5(?—?')/m in dimensionless 
form 



A f 7 2 r 2 + -/V8-7T— |i/>(r)| 2 — 2iV j ^f'^l ^r'] $(r) 



\r — r 
e ip{r) . 



(3) 



In ([3]), e is the chemical potential, and the dimension- 
less quantity 7 denotes the quantum energy of the trap- 
ping frequency in units of the " Rydberg" energy 



7 = hw /E u . 



(4) 



Small values of 7 imply that the effects of the trap- 
ping potential are small compared with the effects of 
the gravity-like interaction, and vice versa for large val- 
ues of 7. In (J3j) we have also assumed N » 1 so that 
the usual prefactor (N — 1) in the Hartree potential can 
be replaced with the total particle number N. Using 
the " order parameter" '5 = ^/N ip instead of the single- 
particle orbital, one can absorb the TV-dependence in (|3|) 
in the wave function 4", and obtains an extended (or, 
for vanishing gravity-like interaction, the familiar) time- 
independent Gross-Pitaevskii equation. 

From (|3|) it would seem that there are three physical 
parameters governing the problem: the trap frequency 
cl>o, given by the dimensionless quantity 7, the particle 
number N and the relative strength a/a u of the scat- 
tering and the gravity-like potential. For the example 
mentioned before one has an estimate of a ~ 10~ 9 m, 
a u ~ 2.5 • 1(T 4 , thus a/a u - KT 6 - 1(T 5 . 
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FIG. 1: (Color online) Phase diagram N versus a/a u for the 
self-binding ground state of monopolar degenerate quantum 
gases. (Explanation see text.) 



However, a central result of the present paper is that 
the physics of degenerate monopolar gases depends only 
on two relevant parameter, viz. 7/TV 2 and N 2 a/a u . To 
see this we note a remarkable scaling property of the 
mean- field Hamiltonian in ((3J: Let tp{r) be a solution of 
the (formal) one-boson problem for a given scaling length 
a/a u and trap frequency 7, 



H mi (N = l,a/a„, 7 )(r) ^(f) = £ ^ 



(5) 



then tp :— TV 3 / 2 ip(r), with f = f/N, solves the iV-boson 
problem for the scaled scattering length N 2 a/a u and the 
scaled trap frequency 7/iV 2 : 



H mf (N, N 2 a/a u , 1 /N 2 ) (r ) ^ (f ) 

with e 



e ip (?) 
N 2 e. 



(6) 



The proof is straightforward and left to the reader. From 
(J6j) follow scaling properties for the mean-field energy, the 
root-mean square radius of the condensate and its peak 
density, respectively: 



E(N, N 2 a/a u ,j/N' 



\(N,N 2 a/a u ,-f/N 2 ) 
Q\(N,N 2 a/a u ),~//N 2 ) 



) = N 3 E(N = l,a/a u ,j) 



= N 4 Q\( N=lta / autl / N 2} 

= N 4 \m\ 2 - (7) 



Iso-surfaces with constant N 2 a/a u and 7/iV 2 form 
planes in the three-dimensional parameter space 
(7, N, a/a u ) on a logarithmic scale. 



B. Self-binding 

In the case of self-binding we are left with one relevant 
parameter, N 2 a/a u . In Fig. [T] we show the phase dia- 
gram N vs. a/a u for self-binding degenerate monopolar 
quantum gases. Note that because of the scaling prop- 
erty, physics is identical on each of the sloping straight 
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lines with N 2 a/a u — constant. Apart from a numerical 
factor, the relevant parameter N 2 a/a v is identical to the 
quantity su used by O'Dell et al. [l2|, but the univer- 
sal nature of this quantity was not emphasized. The two 
asymptotic regimes of self-trapping dubbed " G" (^grav- 
ity") and " TF-G" ("Thomas-Fermi gravity") in are 
characterized by the size of the scaling parameter. For 
N 2 a/a u ^> 1, the kinetic energy is negligible, and self- 
binding results from the balance between repulsive scat- 
tering and gravity- like attraction. For N 2 a/a u <C 1 scat- 
tering is negligible, and self-trapping appears by a bal- 
ance between kinetic energy and gravity-like attraction. 

We note that the "G" regime corresponds to the 
Newton-Schrodinger scheme of quantum mechanics, 
which is a nonlinear variant of quantum mechanics that 
has been investigated in detail in numero us p ublications 
on quantum measurement 111 El 111 111 |M ill 111- 
In that regime, the extended time-independent Gross- 
Pitaevskii equation ^ turns into the Schrodinger equa- 
tion of gravitationally self-interacting quantum parti- 
cles. It is worth noting that interacting monopolar quan- 
tum gases offer an experimental realization of Newton- 
Schrodinger quantum mechanics. 



III. NUMERICAL SOLUTION, RESULTS AND 
DISCUSSION 



We have determined numerically accurate radially 
symmetric solutions of the extended Gross-Pitaevskii 
equation ^ in dependence on the scaling parameter 
N 2 a/a u both for the self-binding case 7/-ZV 2 = and 
for 7/iV 2 7^ 0. To verify the numerical results two dif- 
ferent methods were employed. One was to integrate in 
parallel equation ([3]) and the Poisson equation for the 
gravity-like interaction numerically outward from r = 
by exploiting the initial conditions for the first deriva- 
tives and setting initial values at r = for the wave 
function ipQ and the effective potential produced by the 
gravity-like interaction Vq. The latter was varied via bi- 
section until convergence of the wave function to zero at 
large values of r was attained. The other method was 
an iterative one: the wave functions determined in the 
preceding step are used to calculate the effective poten- 
tial in the next step, and the resulting one-dimensional 
Schrodinger equation is integrated until self-consistency 
is achieved. The iteration is initialized by a reasonable 
guess for the wave function. 

In Fig. [5] we show our results for the wave functions and 
the corresponding self-consistent potentials for the case 
of self-binding for different values of the scaling param- 
eter N 2 a/a u . It can be seen that for increasing N 2 a/a u 
the potentials grow shallower and the wave functions be- 
come more extended. The figure confirms that asymptot- 
icallyall self-binding potentials converge to a 1/r poten- 
tial [2J|. The case of N 2 a/a u = corresponds to the solu- 
tions of the Newton-Schrodinger equation [I^J^IIH H3 • 
As already pointed out by O'Dell et al. jl2|| solutions 
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FIG. 2: (Color online) Numerically accurate self-binding 
ground state s-wave solutions for different values of the scaling 
parameter N 2 a/a u : a) wave functions; b) self-binding poten- 
tials. Both in a) the value of the wave function at the origin 
and in b) the absolute value of the self-binding potential at 
the origin decrease monotonically with the scaling parame- 
ter, from their maximum values at N 2 a/a u — —1.02 to their 
smallest values at N 2 a/a u = 10. Thus, as the scaling param- 
eter grows the binding becomes weaker. In b) the asymptotic 
1/r potential is also shown for comparison. 



also exist for negative scattering lengths, where the con- 
tact interaction, in addition to the gravity-like interac- 
tion, becomes attractive and stability of the condensate 
is established by the equilibrium of the kinetic energy 
of the condensate and the two attractive interactions. 
Fig. [2] shows that for negative scattering lengths the self- 
trapping potentials become ever more binding, until at 
a value of N 2 a/a u w —1.0251 no solutions can be found 
any more, and the condensate becomes unstable with re- 
spect to collapse. This corrects the variational value of 
N 2 a/a u = -3tt/8 « -1.18 given by O'Dell et al. pj. In 
their variational calculation, a Gaussian type orbital was 
assumed, and the mean-field energy of the condensate 
was minimized with respect to the width of the Gaus- 
sian. 

Since we have the numerically accurate solutions at 
hand, we are in a position to check the accuracy of the 
variational results for observables of the condensate ob- 
tained by O'Dell et al. [Hi]. Fig. [3] shows the behavior 
of the total energy of the condensate over seven decades 
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FIG. 3: (Color online) Total energy of the condensate as a 
function of N 2 a/a u : a) on a logarithmic and b) on a linear 
scale. Variational results obtained by minimizing the total 
energy for a Gaussian type orbital yl| are shown by dashed 
lines. 
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FIG. 4: (Color online) Root-mean-square radius of the con- 
densate as a function of N 2 a/a u : a) on a logarithmic and b) 
on a linear scale. Variational results 42J are shown by dashed 
lines. 



of the scaling parameter N 2 a/a u . To cover the range of 
negative scattering lengths, the energy is also given on a 
linear scale in the range around N 2 a/a u ~ 0. The tran- 
sition between the two asymptotic regimes G and TF-G 
around N 2 a/a u ~ 1 is evident from Fig. [3l The com- 
parison with the variational results also plotted in Fig. [3] 
shows that the TF-G regime is well described by the vari- 
ational calculation. It is only in the transition to the G 
regime, and in particular for negative values of N 2 a/a u , 
that sizeable deviations can be observed, up the order of 
10 per cent. 

Observables other than the energy are more sensitive 
to the accuracy of the wave function. We therefore com- 
pare our numerically accurate results with the variational 
results for the root-mean-square radius and the peak den- 
sity of the condensate in Fig. [4] and Fig. O again over 
seven decades of the scaling parameter on a logarithmic 
scale, and on a linear scale around N 2 a/a u w 0. Again 
the transition between the two asymptotic regimes can 
be seen. It can also be recognized that the variational 
results well reproduce the overall behavior of the observ- 
ables. For the extension of the condensate sizeable devi- 
ations occur again in the transition to the G regime and 
for negative values of N 2 a/a u , while for the peak density 
the variational calculation overestimates the correct val- 
ues in the TF-G regime, and underestimates them in the 
G regime, and for negative scattering lengths. Here the 
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FIG. 5: (Color online) Peak density g of the condensate as a 
function of N 2 a/a u : a) on a logarithmic and b) on a linear 
scale. Variational results [l2l ] are shown by dashed lines. 
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FIG. 6: (Color online) a) Bifurcation of the chemical potential 
and b) bifurcation of the total mean-field energy at the critical 
point N 2 a/a u — —1.0251, for self-binding, i. e., vanishing trap 
potential. 



deviations increase up to more than 100 per cent. This is 
understandable since the peak density depends crucially 
on the correct wave function. 



y/N 



FIG. 7: (Color online) a) Bifurcations of the chemical poten- 
tial for nonvanishing values of the trapping potential. The 
case 7 = is shown for comparison, b) Dependence of the 
critical scattering length (the bifurcation point) on the fre- 
quency of the trapping potential. Numerically accurate re- 
sults are given by solid lines, variational results by dashed 
lines. To elucidate the behavior for small values of 7/iV 2 , this 
region is shown in the inset on a logarithmic scale. 



IV. BIFURCATING SOLUTIONS 

A new result of our numerical calculations is that for 
negative scattering lengths there exists a second radi- 
ally symmetric nodeless solution of the extended Gross- 
Pitaevskii equation ((3|). In Fig. [6] the chemical potential 
of the two solutions are plotted as functions of the scal- 
ing parameter N 2 a/a u for 7 = 0. It is evident that the 
critical value of N 2 a/a u — —1.0251 corresponds to a bi- 
furcation point of the eigenvalue spectrum of equation 
©: below the critical point no solution exists, at the 
critical point two solutions appear in a tangent bifurca- 
tion. The total energies of the condensates corresponding 
to the two solutions are also shown in Fig. [5J It can be 
seen that the energy increases from the bifurcation point 
on for the second solution. This would mean that above 
the bifurcation point there exists a collective excited state 
of the condensate in which all atoms occupy one and the 
same nodeless orbital, just like in the true ground state. 

The second solution is in fact present also in the vari- 
ational calculation. It there appears as a second sta- 



tionary (maximum) point of the mean-field energy given 
as a function of the width of the Gaussian type orbital. 
The variational results for the chemical potential and the 
total energy of the second solution are also included in 
Fig. [5] It can be seen that the numerically accurate cal- 
culation is necessary for the quantitative description of 
the bifurcation. 

The second solution even persists in external trapping 
potentials, for any value of 7. The bifurcation diagram 
for two finite values of 7 is shown in Fig. [7J It can be 
seen that with growing 7 the bifurcation point is shifted 
to smaller absolute values of N 2 a/a u . The increase of 
the total energy of the second solution which is evident 
from the figures is a consequence of the fact that the self- 
consistent potentials become more and more binding and 
the wave functions more and more localized which leads 
to a dramatic increase in the kinetic energy. 

What is the physical meaning of the second solution? 

We note, on the one hand, that it corresponds 
to a maximum of the mean-field energy functional. 
Schrodinger's equation, however, and in our case equa- 
tion ^ , follows as the Euler-Lagrange equation of a vari- 
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ational principle which only demands the energy func- 
tional to be an extremum. Thus the fact that the second 
solution corresponds to a maximum of the energy func- 
tional does not preclude it from corresponding to a real 
physical quantum state. On the other hand, the two solu- 
tions are nodeless, and hence nonorthogonal. Obviously 
this is a consequence of the nonlinearity of the extended 
Gross-Pitaevskii equation (J3J): each solution creates its 
own self-consistent potential and thus sees a different 
Hamiltonian. This would seem surprising since the orig- 
inal many-body Hamiltonian is Hcrmitian and linear in 
the wave function, and therefore should possess only or- 
thogonal eigenstates. The nonlinearity of ([3]) is a result 
of the Hartree approximation made for the states. 

In studies of the decay rates in attractive trapped 
Bose-Einstein condensates, with contact interaction only, 
Huepe et al. [28|,[29| have seen similar behavior, i. e., a 
second solution is born in a tangent bifurcation together 
with the ground state. These states also are nonorthog- 
onal. Analyzing the stability of the states Huepe et al. 
have shown that the first excited state out of the two so- 
lutions is unstable with respect to macroscopic quantum 
tunneling. 

This is a strong indication that the second solution 
found in this paper in Bose condensates with gravity- 
like interaction also corresponds to an unstable collec- 
tively excited state. A way to establish this is to lin- 
earize the time-dependent Gross-Pitaevskii equation cor- 
responding to ([3]) around the stationary states and to 
carry out a stability analysis, as was done for the case 
of a pure attractive contact interaction by Huepe et al. 
[Ill [2^| . Alternatively, by choosing a Gaussian ansatz 
with time-dependent widths [3(ij |. equations of motion 
for the widths can be obtained from the time-dependent 
Gross-Pitaevskii equation and analyzed with standard 
stability methods of nonlinear dynamics. Investigations 
along these lines are under way. 

We finally note that there is an analogy with bifurca- 
tions seen in investigations of attractive one-dimensional 
Bose-Einstein condensates on a ring (cf., e. g., [U, HH, 
[33j). There, at a critical value of the ratio of the mean- 
field interaction energy to the kinetic energy, symmetry- 
breaking, soliton-like solutions appear, in addition to 



the symmetry-preserving solution of the Gross-Pitaevskii 
equation, which are lower in energy. By contrast, in the 
example discussed in this paper, both bifurcating solu- 
tions possess the same symmetry. 



V. CONCLUSIONS 

We have re-analyzed Bose condensates with attrac- 
tive 1/r interaction by introducing appropriate atomic 
units which are in particular adapted to the case of self- 
binding. We have thus been able to derive new scal- 
ing properties of such condensates. We have calculated 
numerically accurate results for wave functions and ob- 
servables of self-binding condensates and compared them 
with previous variational results. It turned out that in 
particular at negative scattering lengths the variational 
results become poor and have to be replaced with our 
accurate numerical results. As a novel finding we have 
demonstrated that the critical point where collapse of the 
condensate occurs at negative scattering lengths is in re- 
ality a bifurcation point of the energy functional where 
both the ground state and an excited state merge and 
disappear. We have argued that this second solution in- 
dicates the existence of an unstable collectively excited 
state at negative scattering lengths in degenerate Bose 
condensates with long-range attractive 1/r interaction. 

Critical points, below which collapse of the condensate 
sets in, not only exist in attractive condensates at nega- 
tive scattering lengths [H, [29|, H3, HH, [H, and in the 
monopolar gases with gravity-like interaction discussed 
in this paper but also exist in dipolar gases, in certain 
parameter ranges of the particle number, the scattering 
length and the trap frequencies 0, [H, • Our investi- 
gations suggest that these also correspond to bifurcation 
points. Studies of the bifurcation scenarios in dipolar 
gases are therefore strongly encouraged. 
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